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Abstract. The Navier-Stokes-Voigt (NSV) model of viscoelastic incompress- 
ible fluid has been recently proposed erularization of the 3D Navier- 
Stokes equations for the purpose of direct numerical simulations. In this 
work we investigate its statistical properties by employing phenomenological 
heuristic arguments, in combination with Sabra shell model simulations of the 
analogue of the NSV model. For large values of the regularizing parameter, 
compared to the Kolmogorov length scale, simulations exhibit multiscaling 
inertial range, and the dissipation range displaying low intermittency. These 
facts provide evidence that the NSV regularization may reduce the stiffness 
of direct numerical simulations of turbulent flows, with a small impact on 
the energy containing scales. 
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1 Introduction 

In this work, we study the statistical properties of the three-dimensional 
Navier-Stokes- Voigt (NSV) (sometimes it is written as Navier-Stokes-Voight) 
equations, an incompressible viscoelastic model introduced by Oskolkov in 
[28] , and proposed in pj] and [17] as a smooth regularization of the 3D 
Navier-Stokes equations, for the purpose of Direct Numerical Simulations 
(DNS). 

Throughout the work we consider the NSV model subject to periodic or 
no-slip boundary conditions, and driven by a given force field f . The velocity 
vector field, u(x, t), and the scalar kinematic pressure, p(x, t), are governed 
by the system of equations 

at(u - a^Au) - z/Au + u • Vu + Vp = f, x e fi, 
V ■ u = 0, X G fi, 
u(x, 0) = Uo(x), X G ' 
u(x, t) = X G dQ, or u(x, t) is periodic; 

in the smooth domain Q C M.^, in the case of no-slip Dirichlet boundary 
condition, and with basic periodic domain Q = [0, L]^ C M^, when equipped 
with periodic boundary conditions. Here, a > is a given length scale 
parameter, and u > is a given kinematic viscosity, such that a^/z/ is the 
relaxation time of the viscoelastic fluid. 

It is easy to see from ([1]) that besides having the same steady state solu- 
tions as the Navier-Stokes equations (NSE), the NSV model satisfies formally 
the same infinite time Reynolds averaged equations as those for the NSE, 
suggesting a strong link with the statistical properties of turbulent flows. 

As it was observed above, the NSV model presents an extra length scale 
associated to the viscoelasticity, the parameter a, besides the well known 
Kolmogorov length scale, 1], (see Section 3 for its definition), which is usually 
associated to the smallest scales of motion in turbulent flows. For large values 
of the parameter a, compared to 1], we observe two distinct regions associ- 
ated to the inertial range of the energy spectrum for the NSV model. The 
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first one obeying the celebrated Kolmogorov power law (with anoma- 

lous correction), followed by a second range of length scales, where energy 
condensates, and it is simply equipartitioned. 

The second power-law, however, vanishes as a is decreased, restoring the 
usual Navier-Stokes inertial range regime. This claim is supported by the 
numerical investigation of the regularized Sabra shell model of turbulence, as 
well as by rigorous results reported in the companion paper, [20!, concerning 
the weak convergence of invariant measures of the NSV to strong stationary 
statistical solutions of the Navier-Stokes equations. 

We also present simulations of the Sabra shell model, with the NSV reg- 
ularization term, which is displaying strong damping of dissipation range 
intermittent effects as a is increased, due to a slowdown of the energy trans- 
fer timescales, see Section 3 for more details. We suggest that by tuning this 
parameter, we may attenuate the strong velocity fiuctuations related to the 
intermittent events, reducing thus the the stiffness of DNS of turbulent fiows, 
with only a small effect on the energy containing scales. 

Another main advantage of the NSV model compared to other regular- 
ization and sub-grid scale models used in ocean dynamics, like hyperviscosity 
[23] or a-models [2l[3llll[9], is the fact that in the presence of physical bound- 
aries the NSV model does not require any additional artificial boundary con- 
ditions which cause difficulties and possibly exhibit non-physical behavior in 
applications, such as non-physical boundary layer, see, e.g., 

2 Some facts about the NSV model 

The Navier-Stokes- Voigt (sometimes written as Navier- Stokes- Voight) model 
of viscoelastic incompressible fiuid, ([T]), was introduced by Oskolkov in [25] . 
and pointed out by O. Ladyzhenskaya as one of the reasonable modifications 
of the Navier-Stokes equations, see [22]. In [28], A. P. Oskolkov studied 
and proved its solvability in different functional spaces. We refrain from 
giving technical details about regularity issues here, but we remark that 
these equations behave like a damped hyperbolic system, see [TB]. Therefore, 
solutions do not experience fast (instantaneous) smoothening of the initial 
data, as it is for parabolic systems like the Navier-Stokes equations. 

This fact could prevent the NSV model from being a reasonable modi- 
fication of the Navier-Stokes equations, however, since we are proposing it 
as a model for direct numerical simulations of turbulent fiows in statistical 
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equilibrium, i.e., after the solutions reach the global attractor, we are mainly 
interested in its long time behavior. Indeed, it was proved in [17], that solu- 
tions in the global attractor are smooth, if the forcing field is smooth enough, 
even for initial data satisfying only finite kinetic energy and finite enstrophy 
(i.e. bounded in the Sobolev H^-norm). In particular, in [17], it is shown 
in the periodic case, that if the forcing field is analytic, then the global at- 
tractor consists of analytic functions. This result, in conjunction with results 
proved in [30], proves that if the forcing field is smooth enough, then averaged 
structure functions, with respect to an invariant measure for the NSV flow, 
display exponential decaying tail. One of the main goals of this paper is also 
to present direct numerical simulations, of the Sabra shell model analogue of 
the NSV model, supporting this asymptotic smoothening behavior. 

Let n = [0, L]3. We denote by LP{n), for 1 < p < cx), and H'^iVl) - the 
usual Lebesgue and Sobolev spaces of the periodic functions on fl respec- 
tively. The inner product in the spaces L^(f2), which will be of particular 
interest for us, is given by 

(u, v)= / u(x) ■ v(x) dx, 
Jn 

and its associated norm is defined 

|u| = (u,u)^/l 

Let J-' be the set of all vector trigonometric polynomials on the periodic 
domain Q, and denote 

V = |v9 G : V • ^ = 0, and j ^{x)dx = o|. 

We set H and V to be the closures of V in the and H^{Q) topologies, 

respectively. We equip the spaces H and V with the inner products (■, ■) and 
((■, ■)), and with the corresponding norms |-| and ||-||, respectively. 

The NSV model satisfies the following energy equation for every t e 
[0,oo), 

|(^|u(•,^)|Vy|Vu(•,^)|^) = (f,u(•,^))-HVu(■,^)^ (2) 

Therefore, a positive quadratic conserved quantity in the inviscid, 1^ = 0, 
unforced, f = 0, and periodic or no-slip setting, which was proved rigorously 
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in [5], is 

S^ = l\nf + ^\Vn\\ (3) 
which we call the a-energy. The quantity 

^2 = ^ lu^ (4) 

is the usual kinetic energy, and we remark that it is not conserved for the 
inviscid unforced NSV equations. 

Another conserved quadratic quantity is the a-helicity 

Aq = (u — a^Au, cMr/(u)). (5) 

Because the kinetic energy is not a conserved quantity, the arguments 
used by Kraichnan in [20] (see also [8]) to study the turbulent cascade sce- 
nario cannot be employed directly to the kinetic energy. However, all the 
investigation will be carried out instead to the conserved a-energy, 5*^, and 
conclusions will be further recovered for the kinetic energy, 5*2. This strategy 
was also used in [21 [3l HJ [TOl [131 1211 [25| for studying various a subgrid scale 
models of turbulence. 

We denote by Plh '■ ^ H - the Helmholtz-Leray orthogonal projection 
operator, and by A = —Pi^^A = —A - the Stokes operator subject to the 
periodic boundary conditions with domain D{A) = (if^(fi))^ fl V. The 
operator is a positive definite, self-adjoint, compact operator from H 
into H. Therefore, there exists a complete orthonormal basis of H formed by 
eigenvectors, {'Wj}j>i of A, with associated eigenvalues satisfying < Xj ^ 
oo, when j ^ oo (see, e.g., [TTl [33] for details). 

The term -B(u, v) = Plh((u ■ V)v) is a bilinear form associated with the 
inertial term. Taking the inner product in L^, i.e. in H, of the bilinear form 
with a third variable yields a trilinear form 

6(u,v,w) = (5(u,v),w). 

An important relation for the trilinear form is the orthogonality property 
(see, e.g., [SI [IB [32]) 

6(u,v,v) = 0. (6) 

Let us define the component of a vector field u E H, for a wavenumber 

k, by 

Ufe = PkU := ^ UjWj, 
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where Wj, j > 1, are the eigenvectors of the operator A with the correspond- 
ing eigenvalues A^, j > 1. We also define the component u^'m" by 

Ufc',fc" = Ufc. 

k'<k<k" 

Then, we can write the projected NSV equations 

^(ufc'.fc" + a'^Auk'^k") + i^Auk',k" + B{u, u)k',k" = h',k"- (7) 

Let us now obtain the a-energy budget. Taking the inner product in the 
space H of the equation ([7]) with u^/ j^//, we obtain 



2dt 



6(u, u, Uk'^k") - b{u, u, Ufc/^fc//) + (f , Uk'^k" 
14' -4"] + (f,Ufe/,fc//), 



(8) 



where 

e^(u) = er(u)-erH (9) 
is the net rate of a-energy transfer at k, and 

e^'^(u) = -(5(ufei,fc,Ufei,fc),Ufc_oo) 

represents the net rate of a - energy from the lower modes to the higher 
modes, while that 

e^'^(u) = -(S(Ufc,oo,Ufc,oo),UA:i,fc) 

represents the net rate of a - energy from the higher modes to the lower 
modes. 



3 Averaged Energy Budget 

In this section, we follow [8], [10] and [11] to investigate the energy distri- 
bution scale-by-scale for the 3D NSV equations. Let (■) denote average with 
respect to an invariant measure, fi", for the NSV semigroup (such a measure 
is known to exist for the NSV, see [30]). For the Navier-Stokes equations. 
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assuming that there exists an extensive range of wavenumbers, where the 
viscous dissipation does not play a significant role, one can show that the 
energy simply cascades through these length scales, with rate equals to the 
mean energy dissipation rate for the NSV, = z/(|Vu|^). For the NSV 
equations, a similar scenario holds for the a-energy 

na 1 I |2 , IV7 |2 

•^2 = 2 1^1 + Y l^^l • 
3.1 Energy distribution scale-by-scale 

In the course of investigation of the 5*2 scaling, we will follow methods previ- 
ously used in [21 [3l m [TOl [211 [25] . As it is usual in the studies of homogeneous 
turbulence, we will consider the forcing f with finite number of the eigen- 
modes, i.e. 

f = 5^ f,. (10) 

k<k<k 

Let k" > k' > k. If we take averages in with respect to an invariant 
measure, /i", we obtain the following balance equation 

i^(|Vu,,,.|') = (et'H) - Mn)) + ((f,u,,,„)) = (e^,(u)) - (e^„(u)), (11) 

where e^(u) was defined in (|9]). 

The expression on the right-hand side of the last equality is the mean 
net a-energy transfer in the energy shell [k', k"]. In particular, if we choose 
k" = oo, we obtain 

W|Vu,,^|^) = (e^(u)). (12) 

This expression shows that the net a-energy transfer is positive for every 
k > k. Moreover, if we assume that there exists a range of wavenumbers, 
[k', k"], where the left-hand side of ( fTTl) . i'{\'Vuk',k"f) , is very small, then the 
a-energy transfer is near constant within this range, i.e., 

(e^,(u))~(e^„(u)). (13) 

Defining k^ = {{\Vu\') / {\u\')y/\ we can follow [11] to derive 

- (e^,(u)) - UJ 1^ VA:J ) 
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Indeed, 

^ {eUu)) _ {Ek'<k<k" ^ {k"f{T.k'<k<k" l^fcl') 

" {T.u'<kk'\^k?) \kr) {T.u'<kk'\^k?) 

'{9"(-{£)T' 

Therefore, if k" <^ /c^, then (e^„(u)) ~ (e^,(u)), which means that there is 
no leak of energy in this range. Of course, we cannot expect this condition 
to be fulfilled for every forcing term f . In the last section, where we investi- 
gate numerically the NSV spectrum scenario, we provide Sabra shell model 
simulations satisfying an analogue condition, assuring the cascade scenario. 
We denote by e^, the total energy dissipation rate for the NSV, 

= i/(|VuH. 

Now, we want to investigate the distribution of the inviscid conserved quan- 
tity, the a-energy, scale-by-scale. We define the following characteristic ve- 
locities at scale k: 

uf = (|u.r)v^ 

and 

We denote the characteristic a-energy at scale k by 

sm = (14) 

The characteristic kinetic energy at scale k by 

S2{k) = \{Uff. 
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With this notation, we can write the a-energy as 



(a) 
k ■ 



and the kinetic energy as 

S2 = Y.S2{k) = \Y.{Vl 



{0)n2 



In the inertial range, the a-energy transfer time-scale can be defined as 

^transf^:^_ (15) 

Following arguments used by Kraichnan in [20j, (see also [8], [10] and [3T]). 
in the inertial range, the eddies of size k~^, in average, transfer their charac- 
teristic a-energy to neighboring eddies in the time, t^^^^^^, it takes to travel 
their own length, k~^, i.e., 

,transf_J_ /-.^^n 



where Uk is the characteristic velocity at scale k. Substituting (1161) in (iTSll . 
and setting a = 0, we recover the scaling for the inertial range of the 

S2 structure function, that was theoreticaly predicted for the Navier-Stokes 
equations by Kolmogorov in [T9j (we remark that the k"^/^ scaling for the 
structure function is commonly quoted in terms of its correspondent energy 
spectrum density, which obeys a power law). This sort of argument 

also leads to the double cascading scenario for 2D turbulence described in 
[20| (see also [8]). 

For the NSV case, the situation is complicated by the fact that we have 
two different characteristic velocities, f/^°'' and U^f^. And, in fact, any log- 
convex combinations of them would give us a possible characteristic velocity, 
see, e.g., PSSIEEI. 

Real world turbulent fiows, however, present anomalous scaling, i.e., 
structure functions deviate significantly from the Kolmogorov predictions 
in [19], see, e.g., [T2t [271 132] . This anomalous behavior is present in some 
phenomenological models, as the Sabra shell models introduced in [27]. For 
example, in ^7], the scaling computed for the 5*2 structure function was 0.72, 
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slightly deviating from the Kolmogorov —2/3 scaling. The nature of inertial 
range intermittency is a topic of current intense research in the turbulence 
community, see, e.g., [7] for more details. 

Simulations of the Sabra-NSV shell model, presented in the coming sec- 
tion, clearly display two distinct power-laws for 5^ . For large values of a 
compared to the Kolmogorov dissipation length scale, i] := {v^ j ^a)^^^ ■, we 
observe a range with scaling slightly deviating from the Kolmogorov 
scaling (see Section H] for more details), and another range with a nearly 
power zero scaling, see Figures H] and HI This distribution can be explained 
if we set up the transfer time-scale, 4^^^^^^, in equation (ITBi) . as a function 
of the translational velocity, t/fc ~ '■ 

^transf = ^ ^ {l^a^^fl' 

In fact, substituting the expression above into (ITBi) . we obtain 

Snk)--eTk''/%l + a\Y^^ (18) 

And, therefore, for k -C a~^, we have a range, while that for a ~ 

k~^, we have a power zero range, just as it is observed in the shell model 
simulations in the next section. 

We remark that this scenario of two power laws in the inertial range was 
first proposed in [TU] for the NS-a model, and then for the rest of the a 

models in El la [m Eg. 

3.2 Smallest scales of motion 

The idea of a smallest scale of motion in turbulent flows was introduced by 
Kolmogorov in [19j. This can be obtained by simple dimensional analysis, but 
can also be obtained by comparing the energy transfer time scale, t\f^^^^ and 
the dissipative timescale, i^^^^^P ;= -i^^ leading to the Kolmogorov length 

scale r] = (I'^/ea)^^'^ (see, e.g., [12]). 

In the same spirit, for the NSV case, we define the smallest scale of 
motion, r]^^^ , as the scale where energy transfer time scale, tjj^^^sf^ equals 

the dissipative time scale, t^^^^^^. In order to estimate it, we first obtain 
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explicitly the transfer time scale. Comparing (1141) with fllSp . we obtain the 
following expression for the characteristic velocity at scale k: 

Uf^ ^^J^k'^/^l + a'k^)-^/^ (19) 

Therefore, 

transf 1 il + o'eyl^ , , 

kUf e]I\^l^ 

Thus, when we equate ( 12Ui) with t^^^^^^, we see that for a < t], the smallest 
scale of motion, 77^"^^, is exactly the Kolmogorov length scale, i.e., rj^^'^ = rj. 
However, when a > r/, it is easy to see that 



NSV / \ 1/3 



(21) 



Equation fl2T]) shows that if we choose r/ <^ a, the degrees of freedom of the 
NSV flow are significantly reduced in comparison to the NSE flow, which 
is also a great advantage concerning direct numerical simulations. Figure H] 
supports this fact for Sabra-NSV simulations. 



3.3 Comparison with previous works 

There is no consensus about the correct transfer time scale in previous works. 
In the first studies on the subject, for example in [10], in the context of the 
NS-a model, it was argued that the time scale, t^^^^^^, should be setup by 
the translational velocity, i.e., Ufc ~ U^'''*. This argument was reinforced by 
new large scale simulations by Graham et al. in [15] , although the numerical 
results are mostly inconclusive, due to a still low resolution. 

Recently, Lunasin et al. in [25] and [21], argued that the time-scale should 
be set up by the combination of velocities appearing at the conservation law. 
This leads to the observed power laws for DNS of the 2D Leray-a, in [25] . 
and of the 2D NS-a model, in [21]. 

So far, there is no clear phenomenological explanation for the scaling 
here described for 5*2 . It is in consonance with the arguments for the 3D 
NS-a, obtained in [10], where DNS simulations are still inconclusive, but 
in dissonance with the arguments of the 2D simulations in [211 12S]5 where 
high resolution DNS were employed, and clear power laws were presented. 
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It is reasonable that the form of the nonhnear term, responsible for the 
energy transfer mechanism, should play a role in the transfer time scale of 
the conserved quantity. However, none of the former arguments contemplate 
it, and a robust explanation for the observed scaling is still lacking. 

We also observe that in the shell model simulations in the next section, 
the distribution of kinetic energy scale by scale, 5*2, presents an inertial range 
with three distinct power laws. Indeed, the log-log plot of 5*2 in Figure H] 
follows the shape of S"^, with its two power laws, until it reaches wavenumbers 
obeying a third power law, with scaling k~^, just before the dissipation range 
takes place. We remark that this intermediate range, with a power zero 
spectrum, which is reminiscent from the a-energy spectrum, has never been 
theoretically predicted or numerically computed in any of the former works 
about a models, see, e.g., [21 El I2 [lOl [161 [Ml [25] . 



4 Sabra-NSV shell model simulations 

Due to the extraordinary complexity of the hydrodynamic equations, many 
simplified models, based on the phenomenological theories of turbulence, 
have been proposed in order to investigate the statistical scenario of complex 
flows. In this section, we use a modification of the Sabra shell model of 
turbulence, which describes the evolution of complex Fourier-like components 
of a scalar velocity field denoted by Un G C. The associated one-dimensional 
wavenumbers are denoted by where the discrete index n is referred to 
as the "shell index". The equations of motion of the Sabra shell model of 
turbulence were introduced in [27], and they have the following form 



= i{akn+lUn+2Un+l + bknUn+lU*„_i - ckn-lUn-lUn-2) " l^klun + fn, (22) 

for n = 1, 2, 3, . . . , with the boundary conditions U-i = uq = 0. The wave 
numbers kn are taken to be 

kn = koX"", (23) 

with A > 1 being the shell spacing parameter, and ko > 0. Although the 
equation does not capture any geometry, the scale L = k^^ is frequently 
considered as a fixed typical length scale of the model. In an analogy to the 
Navier-Stokes equations z/ > represents a kinematic viscosity and /„ are 
the Fourier components of the forcing. 
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Figure 1: Log-log plot of Sabra-NSV simulation with parameters v = 10^^, 
and a = 10~^. (+) 5*2 (A;) - Characteristic a-Energy at scale k. (o) 5*2 (/c) - 
Characteristic Kinetic Energy at scale k. 



We consider the following regularization of the model associated to the 
NSV equations, which we denominate as the Sabra-NSV shell model 

~ir = 1 , \, 2 iakn+lUn+2U*n+l + hk^Un+lKi-l - ckn-lUn-lUn-2) 

dt 1 + a^k^ ^ ^2^^ 

^^n , fn 

:Un + 



l + a^kl l + a^kl'' 



for n = 1, 2, 3, ... , for fixed a > 0. 



The three parameters of the model a, h and c are real. Following |27j, we 
require that in the inviscid (z/ = 0) and unforced (/„ = 0, for all n) case, the 
model should have at least one formal quadratic positive definite quantity 
to be invariant. Such a quantity will represent the a-energy in the system. 
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Figure 2: Log-log plot of the Sabra-NSV simulation with parameters 
10~^, and a = 10~^. (+) 5*2 (fc) - Characteristic a-Energy at scale k. 
5*2 (/c) - Characteristic Kinetic Energy at scale k. 
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Indeed, in order to require that the a-energy 

oo 

E = J2i^+a'kl)\ur,\', (25) 

n=l 

will be formally conserved, we assume the following relation between the 
parameters of the model, which we will refer as an "energy conservation 
assumption" 

a + 6 + c = 0. (26) 

Moreover, in the inviscid and unforced case the model possesses another 
formal quadratic invariant 

w=E(^) (i+«^aKr- (27) 

n=l ^ ^ 

For f < this quantity is not sign definite and thus it is common to associate 
it with the "helicity" - in an analogy to the three-dimensional turbulence. 

Remark. For a = 0, that is, when we are in the pure Sabra regime, the 
above mentioned conservation laws are only formal, due to a possible lack 
of regularity of solutions of the inviscid Sabra shell model. However, by 
following [3] (see also [6]), we can prove, when a > 0, global existence and 
uniqueness for the inviscid Sabra-NSV shell model, and therefore, it can be 
proved that E and W are rigorously conserved in the inviscid and unforced 
case. 

Without lost of generality we may assume that ko = 1. Next, by rescaling 
the time 

t at, 

and using the "energy conservation assumption" (126|) we may set 

a = 1, b=-9, c = e-l. (28) 

Therefore, the Sabra shell model is in fact a three-parameter family of equa- 
tions with parameters u > 0, 6, and A. We are interested in the case where 
the shell sizes grow geometrically (see (123|1 ). therefore we limit ourselves to 
A > 1. 
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The three-dimensional parameters regime corresponds to < ^ < 1, as 
W is not sign definite (see, e.g., [H], [SZ])- In that regime we can rewrite 
relation fl27|) in the form 

oo 

W = 5^(-irA:f(l + a^A;^)K|2, (29) 

n=l 

for 

(3 = -\og,\e-l\. (30) 

In our simulations, we set 6* = 1/2, A = 2, fco = 1, and with constant 
forcing, of order one, in both real and imaginary parts of the second and 
third modes. 

We first comment on the results derived in the last section for the charac- 
teristic energy distribution scale by scale. In shell models, we refer to these 
quantities as the second order structure function, and define them by 

and 

S!^{K) = {\unf + a'e\unf). 

Figure 1 and Figure 2 show the results of simulations for u = 10~^, and for 
a = 10~^, and a = 10~^, respectively. They show two clear distinct power- 
laws, one with a scaling of and another one with a constant scaling. 
Figure 3 shows the effect of varying the parameter a for fixed u = 10~^. As 
expected, the secondary power law becomes less prominent as a is decreased. 

Now, we comment on some issues about intermittency. One of the main 
characteristics of turbulent fluid flows, that is present in the Sabra shell 
model, is its dissipation-range intermittency. This is characterized by violent 
fluctuations of very short duration in the energy dissipation rate, 

= \'^n\ )■ 

n 

Small velocity fluctuations in high wavenumbers play a key role in this phe- 
nomenon, see, e.g. [12], and the rise of such fluctuations imposes severe 
challenges to the direct numerical simulations of turbulent flows. Therefore, 
the attenuation of its effects flnds many applications. 
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Figure 3: Log- log plot of the characteristic kinetic energy scale-by-scale for 
the Sabra-NSV simulations with v = 10~^, and different values of the pa- 
rameter a. 

Looking carefully at the equation (IMl) . we observe that as kn becomes 
large, the term 1 -|- a^fc^ in the denominator will damp the velocity fluc- 
tuations. Therefore, the energy dissipative intermittency must be signifi- 
cantly attenuated for large values of the relaxation time parameter. Fig- 
ures 4, 5 and 6 present the energy dissipation fluctuation signal, e'^/e^ = 
{^n>i^l\'^ri{t)f)/{Y,n>i^l\'"'nit)f), ^OT different settings. The intermit- 
tency becomes strongly attenuated as we increase the length parameter a. 
This fact might find real world applications, and worth further investigations. 

Another main characteristic of the Sabra model, as observed in [27], is 
that it exhibits an inertial range, with moments of the velocity depending 
on kn as power laws with nontrivial exponents, i.e., {\un\'^) ~ kn^'' , where 

depends nonlinearly on q. From now on, we refer to these moments as 
structure functions. For even q = 2m, we use the usual definition 



S2m{kn) = (|U 



2m 



'n 



For odd q = 2m -|- 1, we use the following definition 



S2m+l{kn) = Im{Un-^iU„Un+l \u 



2m-l 



)■ 
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For small values of the length parameter a, i.e., for a rj, where rj is the 
Kolmogorov length scale, we do not observe any significant deviations for 
the exponents of the structure functions from the pure Sabra model, with 
a = 0. In Figure 7, we show some structure functions for a simulation with 
u = 10^^, and a = 10"''. The values of the anomalous exponents are listed 
in the tigure. 

For larger values of the length parameter a compared to the analysis 
of the anomalous exponents is complicated by the presence of the double 
power law, and since we are concerned with the small relaxation time regu- 
larization, we will not present any detailed analysis of it here. 

However, we want to announce an interesting fact that we plan to in- 
vestigate in a forthcoming work. Self-similarity is a crucial hypothesis in 
K41 theory, see, e.g., [12]. Intermitent dynamics in turbulence, however, 
is inconsistent with the self-similarity hypothesis, leading to several modi- 
fications of this hypothesis, see, e.g., |12], for more details. One can mea- 
sure the departure from self-similarity in the inertial range by looking at 
the flatness of the signal in this range, see [12]. The flatness is defined as 
F{k) = {Si{k)) / {S2{k)'^) , and a departure from a constant value means a 
lack of self-similarity. For turbulent-fiuid fiows, like in shell models, this 
term deviates significantly from a constant as k = kn increases, even within 
the inertial range. 

The nature of inertial range intermittency in shell models is not well un- 
derstood. The infiuence of high wavenumbers strong velocity fiuctuations 
might play a role, just like in real turbulent fiows, as suggested by Landau 
in his famous criticisms of K41 theory, see, e.g., [12]. Because for the NSV 
model, we have strong damping of high wavenumbers velocity fiuctuations, 
we calculated the fiatness of the Sabra NSV model to investigate the rela- 
tionship between the two kinds of intermittency. 

We observe that the inertial range intermittency is significantly reduced as 
the parameter a becomes large when compared to the Kolmogorov dissipation 
length scale, i]. For example, for the pure Sabra model with u = 10"^, we 
observe a fiatness of the order F{kn) ~ while that for the Sabra-NSV 
with same viscosity, and values of a = 10~^, and a = 10~^, the fiatness is, 
respectively, F{kn) ~ k^°^, and ~ 

For u = 10~^, the departure was from F{kn) ~ k^^^, for the pure Sabra 
model, to F{kn) ~ for a = 10-^ and to ~ A;0 °02 for a = 10"^ 

Another interesting remark is that we fit the first power law, and we observed 
a trend towards the KAl scaling with less and less anomaly, as the relaxation 
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Figure 4: Log-log plot of the energy dissipation fluctuation e^/e^. Parame- 
ters: V = 10~^. a = 0. 



time was increased. 

Genuine fluctuations for odd order structure functions, see [27] , prevented 
us from performing more detailed analysis of this pattern, and we are cur- 
rently investigating it by using a special stochastic forcing that eliminates 
this complication. We expect to report the results of this investigation soon. 
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Figure 5: Log-log plot of the energy dissipation fluctuation e^/e^. Parame- 
ters: V = 10-^ a = 10-^. 
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Figure 6: Log-log plot of the energy dissipation fluctuation e'^/ca- Parame- 
ters: u = 10"^. a = 10"^ 
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Figure 7: Log-log plot of the structure functions, S2, S3, S4 and Sq for a 
simulation with u = 10~^, and a = 10~^. One can also find the corresponding 
scaling of the structure functions in the inertial range. 
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